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SUSCEPTIBILITY CORRECTION AND CONTRAST ENCODING 
WITH SIMULTANEOUSLY REFOCUSED MR IMAGES 

10 Field 

This patent specification is in the field of Magnetic Resonance Imaging 
(MRI) and more particularly relates to improving imaging by removing or reducing 
imaging artifacts such as those caused by susceptibility effects and to other 
15 improvements in MRI. 

Background 

Susceptibility artifacts in MRI images are believed to be due to local 
20 inhomogeneities in the main magnetic field, and to be caused by local 
anatomical or physiological features. For example, undesired local 
magnetic field gradients can occur in the head in regions near air-tissue 
interfaces near the frontal sinuses and can cause MRI signal loss in the 
orbital gyrus and in the inferior temporal gyrus and inferior frontal cortex. As 
25 pointed out in Yiping Du, Journal Review: Z-Shimming in MRI, 1/14/00, 

incorporated herein by reference, some of the susceptibility effects may be 
corrected by post-processing but others are believed to be unrecoverable. 
Susceptibility artifacts and corrections are discussed in Constable, RT. and 
Spenser, D.D., Composite Image Formation in z-Shimmed Functional MR 
30 Imaging, Magnetic Resonance in Medicine 42:1 10-1 17 (1999), and in 
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Glover, G.H., 3D z-Shim Method for Reduction of Susceptibility Effects in 
BOLD MRI, both of which are hereby incorporated by reference. 

Despite the progress discussed in the material cited above, and in 
the material cited in the Journal Review by Yiping Du, it is believed that a 
5 need still exists for a more effective reduction or elimination of susceptibility 
effects. 

An echo planar pulse sequence (EPI) is known in MRI and involves 
using a slice- selective RF excitation pulse applied with a slice-selective 
gradient to define an MRI signal within a slice volume. The signal is 
10 refocused by multiple (n) switchings of the read gradient, during which time 
n signals are recorded — one signal refocused after each read gradient 
switching. 

A Burst sequence is also known, where a multitude of hard RF 
pulses are beamed into the patient during the application of a defocusing 

15 gradient. The signals are refocused by a single refocusing gradient, during 
which time the signals are also phase-encoded for a second axis of spatial 
resolution. The time of applying plural hard RF pulses, which are not slice 
selective, is much shorter than the time of gradient switching in EPI, so that 
the Burst method is faster than EPI. However, the SNR of the Burst method 

20 is much lower than of the EPI method, as the signal of the Burst sequence 
is reduced by fragmentation of the total available signal magnetization 
proportional to the number of hard RF pulses. 

A hybrid of Burst and EPI also has been proposed - see Heid U.S. 
Patent No. 6,034,528, incorporated herein be reference. After a Burst 

25 sequence, a read gradient is switched in polarity and the echoes in the train 
created by the hard RF pulses are again refocused. The signals generated 
by the hybrid sequence are used to generate a 2D FT image by using RF 
saturation pulses outside of the imaged slice volume. The Burst method 
uses hard RF pulses that do not have slice selectivity by imposing a 

30 frequency modulation on the RF pulse during the application of a slice 
selective (Gss) gradient with an opposite polarity pulse on the same Gss 
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axis. The hard RF pulses are applied in much shorter time than the slice 
selective RF excitations, as required to generate the Burst signals very 
quickly to avoid T2 decay of signal. The gradient applied during the hard 
RF pulses of a Burst sequence are on the read gradient axis (Gr) and are 
5 used to defocus the magnetization of each RF pulse differently so that they 
refocus at different times on the switched polarity Gr. Further different from 
the below described z-shim-SIR sequence, the Burst sequence does not 
apply gradient pulses between the hard RF pulses, only a single continuous 
duration of a single Gr gradient that is on at all times during and between 

10 the RF pulses. Further different from the z-shim-SIR described below, the 
Burst sequence combines all of the signals generated by the sequence to 
make one 2D image, or to make a 3D image by either 2D FT or 3D FT. 
Also different in making a 2D image, the Burst sequence saturates signals 
in two regions on either side of a desired image plane; since slice selectivity 

15 is not achieved by using the hard RF pulses. 

Also known is a Simultaneous Image Refocusing (SIR) pulse 
sequence, described in co-pending patent application Ser. No. 09/477,783 
filed on December 31, 1999" in the name of David Feinberg, a co-inventor in 
this application, and also described in DA Feinberg, TG Reese, VJ Wedeea, 

20 ISMRM 2000, p. 681. The 1999 application and the 2000 publication are 
hereby incorporated by reference in this patent specification. In the SIR 
pulse sequence, each of a number of slice volumes is given a respective 
different phase history. The process then simultaneously refocuses these 
slice volumes, and acquires the MRI signals from them at respective 

25 different times related to their respective different phase histories. For 
example, two consecutive 90E RF excitation pulses at respective 
frequencies create magnetization in two slice volumes s1 and s2 of a body, 
such as by using the known MRI practice of employing a main magnetic 
field and a slice selection gradient field acting on a body. A first dephasing 

30 pulse, on a readout axis Gr, is applied after the first but before the second 
90E RF pulse, thereby encoding the magnetization of the first slice volume 
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(s1) but not of slice volume s2 since the second slice volume magnetization 
has not yet been created. A second dephasing pulse, on the same axis Gr, 
is applied after both slice volumes have been magnetized by their 
respective 90E FR pulses and, therefore, encodes the magnetization of both 
5 slice volumes equally. This gives the two slice volumes respective different 
phase histories relative to the Gr direction. If the subsequent readout 
pulses on the Gr gradient axis are of unit area, and the first and second 
dephasing pulses on the Gr axis have areas of 0.50 and 0.25, respectively, 
the magnetization pathways they produce in the two slice volumes will 

10 refocus echoes (MRI signals) within opposite halves of each read period, 
alternating with each EPI (echo planar imaging) gradient switching, due to 
the respective different phase histories of the two slice volumes relative to 
the read gradient direction. Thus, in one example an SIR EPI sequence 
applies two RF excitation pulses on two different slice planes using different 

1 5 frequency offsets with slice-selective gradients. Application of gradient 
pulses on the read gradient axis, after each RF pulse, causes a different 
amount of defocusing of the signal created by the two RF pulses. The 
application of an EPI type of multiple switched read gradients causes 
refocusing of the two different slice signals to occur at two different times on 

20 each switched read gradient. Time reversal of signals further separates the 
signals from the two different slice planes onto two different areas of the 
acquired k-space data, which can be digitally separated into two smaller k- 
space data sets of separated data from each slice. Standard 2D FT 
(Fourier Transform) gives two slice different image for positions obtained 

25 from the single SIR EPI echo train. 

While EPI and SIR processes have made significant contributions to 
MRI, it is still desirable to provide improved MRI processes, and various 
features of the new processes described below are directed to that end. 

30 Summary 
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One object is to extend correcting for susceptibility signal loss caused 
by susceptibility defocusing. To this end, a new process creates two or 
more images from MR signals for the same slice volume, which images 
have different characteristics, and combines them into a single image that is 
5 susceptibility-corrected in a way superior to the known approaches. 

In an example of 2D FT EPI images, the new MRI process acquires 
MRI signals (spins) available in a planar volume during a T2 relaxation 
period. It modifies the known EPI sequence to refocus signals created from 
plural RF excitation pulses that are applied with slice-selective gradients. 

1 0 The multiple magnetizations of spins (m) created by m number of RF pulses 
are each encoded with a different area of gradient pulse to compensate for 
a different defocusing effect of the susceptibility-created local gradients that 
are time invariant and dependent on the object structure and composition. 
The differently gradient-encoded m magnetizations are then refocused by n 

1 5 EPI read gradient switchings and n phase-encoding gradient pulses to 
encode 2D FT images. The pulse sequence generates m 2D FT n x n 
matrix k-space data sets, with each of the corresponding images having a 
different gradient encoding to correct for susceptibility defocusing of MRI 
signals. One of these images has no gradient correction and creates an 

20 image that has regions of signal loss due to the susceptibility. The other m-1 
number of images have MRI signals in the region of susceptibility, according 
to the refocusing of the defocusing by susceptibility by an equivalent but 
opposite sign phase shift on spins created by the amplitude-time area of the 
gradient pulse. As with the known EPI technique, a time reversal is 

25 performed on each sequential group of m signals after every other read 
gradient switching or after alternating read gradient switchings. The m 
images are created of the same image volume at essentially the same time, 
essentially simultaneously from the same switched read gradient sequence. 
Therefore, they experience essentially the same physiologic changes in 

30 blood flow and due to any temporally varying physiologic function including 
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respiratory motions and cardiac pulsations and neurologic stimulations as 
used for functional MRI studies with BOLD contrast. 

The m images are combined using any suitable digital method, and 
most simply by addition after magnitude display of the data. The combined 
5 image created from the addition or other combination of m images differs 
from a conventional EPI image in that the areas of signal loss caused by 
susceptibility are filled in with signal recovered in the images using 
conditioning gradient pulses. The new pulse sequence creates an image 
with reduced or eliminated areas of signal loss from susceptibility, areas of 

10 signal loss that are present in known EPI images. The new pulse sequence 
thus shows a marked improvement over conventional EPI images by 
recovering image signal in the regions of the temporal lobes and frontal 
lobes of the brain, where the susceptibility errors from air in frontal and 
mastoid air sinuses typically appear. 

1 5 Various aspects and variants of the new process are applicable to 

sequences other than the known EPI sequence and to achieving results 
other than suppression of susceptibility artifacts, as discussed in connection 
with examples described in detail below. 

20 

Brief Description of the Drawing 

Fig. 1 illustrates a known EPI pulse sequence using RF excitation " 
and generating one MRI signal per switched read gradient. 
25 Fig. 2 illustrates a known SIR EPI pulse sequence using multiple 

slice-selective RF pulses, where Gss is applied during each RF pulse that 
has modulations, and the RF pulses create magnetization on different slice 
volumes by offsetting the frequency of the RF pulse. 

Fig. 3 illustrates a new, z-shim-SIR EPI pulse sequence where the 
30 same slice volume is irradiated by the slice-selective RF pulses, and 

conditioning gradient pulses are applied between the RF pulses. The Gc 
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conditioning gradient pulses can be applied on each gradient axis, Gr, Gp or 
Gs. 

Fig. 4 illustrate a result from using a sequence such as in Fig. 3, 
where there are two intermediate MR images of the same slice volume but 
5 one is affected by the conditioning pulses while the other is not, and the two 
intermediate images are combined into a single image with reduced or 
eliminated susceptibility dropout. 

Fig. 5 illustrates an example of timing of events in a z-shim-SIR EPI 
process. 

10 Fig. 6 illustrates k-space data collected from the pulse sequence of 

Fig. 3, and a separation of k-space data and the use of FT to create 
differently conditioned images of the same slice volume. 

Fig. 7 illustrates simultaneous acquisition of multiple 2D phase- 
encoded lines with SIR EPI. 
15 Fig. 8 illustrates simultaneous acquisition of multiple 3D phase- 

encoded partitions with SIR EPI. 

Fig. 9 illustrates echo-spacing in k-space. 
Fig. 10 illustrates phase-cycled z-shim SIR EPI. 
Fig. 1 1 illustrates a pulse sequence for a 3-image process from two 
20 RF pulses. 

Fig. 1 2 shows MR images from the sequence of Fig. 1 1 . 
Fig. 13 illustrates a pulse sequence involving velocity encoded 
pulses. 

Fig. 14 illustrates a pulse sequence involving multiple correction 
25 pulses. 

Detailed Description 

Before describing examples of the new, z-shim-SIR EPI and other 
30 process, it may be helpful to describe known EPI and the SIR processes to 
set a background. Referring to Fig. 1, an EPI process involves applying an 
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RF pulse a° concurrently with a slice-selective gradient pulse on a Gs axis 
that is followed by a half-width slice selective refocusing gradient of opposite 
polarity. This is followed in time by readout gradient pulses on a Gr axis 
that alternate in polarity, starting with a half-width pulse. The half-with pulse 
5 on the Gr axis coincides in time with a pre-phasing negative pulse on a Gp 
axis. Blipped phase-encoding pulses coincide with the start and end of 
each negative pulse on the Gr axis after the half-width pulse on Gr. The MR 
echo signals appear after polarity switchings of the Gr pulses, and typically 
are used to reconstruct and display or otherwise use a conventional MR 
10 image. 

The SIR EPI process illustrated in Fig. 2 produces MR echo signals 
for two slice volumes, s1 and s2, in a single echo train. Two 90° RF pulses 
are applied, F1 and F2. A first slice-selective gradient is applied on the Gs 
axis at the same time as F1, followed in time by a half-width slice-selective 

15 refocusing pulse of opposite polarity on Gs. The second RF pulse, F2, 

coincides in time with a second slice-selective gradient pulse on Gs, flanked 
by two half-width slice-selective refocusing pulses of opposite polarity on 
Gs. Thus, two different slice volumes are defined. A half-width, negative- 
polarity, conditioning gradient pulse on the Gr axis is halfway between the 

20 two RF pulses and the Gs pulses, and a quarter-width, negative-polarity, 

conditioning gradient pulse is on the Gr axis. The quarter-width pulse on Gr 
coincides with a pre-phasing pulse on the Gp axis. This gives respective 
different phase or conditioning histories to the two slice volumes, as the 
half-width gradient pulse on Gr affects only one of the slice volumes but the 

25 quarter-width conditioning pulse affects both slice volumes. Then, full-width 
read gradient pulses of alternating polarity follow on the Gr axis, and a 
blipping phase-encoding pulse on the Gp axis matches each polarity 
transition of the full-width gradient pulses on the Gr axis. Two MR echo 
signals appear during each of the full-width read gradient pulses, one for 

30 each slice volume, thus achieving a number of benefits over conventional 
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EPI, such as faster acquisition of MR signals for reconstructing two slice 
images. 

One example of the new process is illustrated in Fig. 3 and extends 
SIR EPI into a z-shim-SIR EPI process by adding conditioning gradient 
5 pulses between two RF pulses labeled 60° and 90°: (a) a gradient pulse 
GsC on the Gs axis; (b) a gradient pulse GrC on the Gr axis; and (c) a 
gradient pulse GpC on the Gp axis. Unlike the SIR example earlier 
discussed, the two slice-selective gradient pulses that coincide in time with 
the two RF pulses, Gs1 and Gs2, select the same slice. Also unlike the SIR 

10 example, there are two pre-phasing gradient pulses on the Gp axis, a 

positive pulse GpP1 that coincides in time with the negative-polarity half- 
width pulse GsP2 on the Gs axis, and a negative-polarity GpP2 pulse that 
follows in time the quarter-width pulse GrP2 on the Gr axis. The other 
gradient pulses on the Gr and Gp axis match those in the SIR example. In 

15 particular, Gs1 is a slice-selective pulse for the slice volume, as is Gs2. 
GsR1 and GsR2 are slice-selective refocusing gradients. GrP1 is a pre- 
phasing pulse for Gr to shift the position of signals s1. GrP2 is a pre- 
phasing pulse affecting the positions of signals s1 and s2 on the readout 
period. Gr+ and Gr- are positive and negative polarity readout gradient 

20 pulses. GpP1 and GpP2 are pre-phasing gradient pulses on the phase- 
encoding axis, and GpBlip are blipped phase-encoding gradient pulses. 

In Fig. 3, two sets of MR echo signals are obtained, a set labeled s1 
and a set labeled S2. Unlike the SIR example, here the two sets are for the 
same slice volume. However, they contain different information about that 

25 slice volume. The conditioning gradient pulses affect one of the two sets of 
MR signals but not the other. Thus, two MR images can be reconstructed 
for the same slice volume, using for example conventional FT processing for 
each. One image is from the set of MR signals affected by the conditioning 
pulses and one from the set that is not so affected. Stated differently, one 

30 of the MR images is corrected, but the other one is not. When the two MR 
images of the same slice volume are combined, by adding them or using a 
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more sophisticated combination of the two, the resulting combined MR 
image is superior - for example it suppresses susceptibility dropout. 

Fig. 4 illustrates an example of such MR signal acquisition and the 
resulting two initial MR images that are for the same slice volume (but 
5 contain different information), and the MR image resulting from combining 
the two images, in this case by simple addition. The left image is the image 
that is not affected by the conditioning gradient pulses; the center image is 
the one affected by those pulses, and the right image is the result of 
combining the left and center images and demonstrates a reduction or 

1 0 elimination of susceptibility dependent signal dropout. 

Fig. 5 illustrates timing of events in the process of Fig. 3. The EPI 
sequence is modified by applying at least one additional RF excitation pulse 
prior to the switched gradient readout of signals. Each RF excitation pulse^ 
creates signal from the same slice volume. A conditioning gradient is 

1 5 applied on the read gradient axis during the period between the RF 

excitations to encode a different signal defocusing on signals from each 
respective RF excitation. The amount of signal defocusing is determined by 
the polarity and gradient pulse area (A), and the amplitude-time relationship 
(y- Gr- Tp) where y is the gyromagnetic constant, Gr is the amplitude of the 

20 gradient pulses applied on the gradient read axis (milliTesla/meter) and Tp 
is the duration (millisecond) of the gradient pulse. After the time of the RF 
pulses and the defocusing gradient pulses, a multitude of switched readout 
gradients are applied to create multiple additional signal refocusings as in 
an EPI echo train sequence. The switched readout gradients refocus the 

25 signal from the RF pulses at different times during the readout gradient. The 
time of each signal refocusing is determined by the net sum of the 
defocusing gradient experienced by the signal. The MR signal created by 
the first RF excitation pulse experiences the net sum of all defocusing 
gradient pulses applied after the first RF excitation pulse. The MR signal 

30 created by the second RF excitation pulse does not experience the 

defocusing gradient pulses applied prior to the second RF excitation pulse. 
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Therefore the first and second MR signals differ in defocusing by the 
conditioning gradient pulse applied between the respective first and second 
RF excitation pulse. This difference in defocusing of the MR signals causes 
the signals to be refocused at different times or locations on the read 
5 gradient pulse. The MR signal is refocused at a time when the time integral 
of the readout gradient creates a phase refocusing equal to the net 
defocusing of the signal. After each read gradient period, there is a 
switching of the polarity of the read gradient, similar to known EPI 
sequences, which continues the process of refocusing additional signals on 

10 each subsequent read gradient. The timing of the first and second signals is 
reversed in each switched read period. At the end of the EPI echo train 
sequence, the signal from each RF excitation is interleaved into every other 
signal of the collected data. After time reversal of the signals from every 
other read period, the two sets of signals from different RP excitations are 

15 on opposite sides of k-space on the kr axis. 

As illustrated in Fig. 6, the k-spaces of two separate images are 
obtained by dividing the data matrix around the center of the Kr-axis. The 
resulting two k-space sets of signals are simultaneously phase encoded in 
the echo train by gradient pulses (Gp) for a second axis of resolution either 

20 using a blipped or continuous Gp gradient waveform during the echo train. 
Standard 2D Fourier Transform (FT) of each of the k-space data sets is 
made to produce the two images at the identical slice position, or some 
other suitable image reconstruction method can be used. 

Now that it is explained how two different images of the same slice 

25 volume are created from two different RF excitations using one sequence of 
multiple switched Gr gradients, as in an EPI sequence, it can be further 
understood how these two images can be encoded separately to correct for 
susceptibility. Any suitable type of gradient or RF pulse used for a 
susceptibility correction scheme that is applied after the first RF pulse will 

30 encode the signals of the first RF pulse. However, just as the signal created 
by the second RF excitation pulse does not experience the defocusing 
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caused by Gr2 gradient pulses applied prior to the second RF excitation 
pulse, the second signal does not experience the correction or conditioning 
scheme pulses. Similarly, a net sum of correction or conditioning scheme 
pulses is experienced by the signal beginning at a time in the sequence 
5 when the respective RF excitation pulse creates the signal, not' including 
sequence time prior to the RF pulse. The term "phase history" is introduced 
to define this net accumulated defocusing of the gradient pulses applied 
during the time of the signal, which excludes pulse sequence effects prior to 
the RF excitation pulses that create the signal. 

10 Susceptibility correction schemes involve the application of gradient 

pulses to compensate or refocus the linear components of the local 
gradients (G L ) created by susceptibility. Signal loss is caused by 
susceptibility. This effect is due to local magnetic field gradients in the head 
in regions near air - tissue interfaces. These regions occur near the frontal 

15 sinuses and cause signal loss in the orbital gyrus, and the inferior temporal 
gyrus and inferior frontal cortex. The defocusing caused by the local 
gradients of susceptibility is independent of other gradient pulses applied 
during the pulse sequences, however the net defocusing due to 
susceptibility is largely determined by the TE of the sequence, the T2*, as 

20 (G L ■• TE) = D* where the defocusing of the signal (D*) is the contribution 
made by linear components of the local gradients of susceptibility. D* is 
refocused (R*) by the combined effect of specific compensation gradient 
pulses (Gc) of duration Tc, where R* = Gc Tc, such that there is refocusing 
of the signal loss by D* by the relationship (D* + R* = 0). In the same 

25 image that has this susceptibility correction by Gc gradient pulse, the signal 
in other regions of the slice that do not have local susceptibility effects is 
defocused by Gc. These regions therefore experience a loss of signal. 
Therefore, the new process creates at least two images, one which has Gc 
and R* correction of susceptibility and a second image that does not have 

30 R* correction. A combination of the two images creates a net image in 

which there is recovery of signal in regions where the defocusing effect of 
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local susceptibility gradient is nulled or refocused by the applied Gc gradient 
pulse. The net effective Gc gradient pulse is applied between the two RF 
excitation pulses, so in effect the image created by the first RF excitation 
pulse has the recovery of signal in regions of susceptibility local gradients. 
5 The image derived from the second RF pulse does not experience the Gc 
pulse, and therefore it has a signal loss in regions of susceptibility but 
normal signal in all regions of the slice where there is no susceptibility 
effect. In a single shot echo train, these two images are acquired. These 
two images are combined to create a final image that has signal recovery in 

10 regions of high susceptibility where the rephrasing of Gc pulse equals the 
defocusing of the local susceptibility gradients. At regions where D* does 
not equal R* f there is signal loss. Nevertheless, by appropriate adjustment 
of Gc, much of the signal loss in EPI images is recovered in the corrected 
image from the combination of two simultaneously acquired data sets. 

15 A variant of the above pulse sequence is to encode 3 (or more) 

images in each sequence. The method can then be extended so that two (or 
more) images with different R* values are encoded and one image has no 
R*. One method of creating 3 different signals is by using 3 RF excitation 
pulses. Another way to create a third image is to create a spin echo signal 

20 by refocusing the signal with two RF pulses. Some of the signal created 
from the first RF pulse is refocused by the second RF pulse to create a 
second set of signals with a different phase than the remaining signals 
created by the first RF pulse. This can be extended to more than three 
images. 

25 A different variant of the above described method and technique is to 

correct for susceptibility and main magnetic field inhomdgeneity errors in 
MR images that show velocity of motion of spins. An example of a suitable 
pulse sequence is illustrated in Fig. 13. The multiple images of the same 
slice volume are encoded with different velocity-dependent phase encoding 

30 pulses. By combining (e.g. adding) the velocity phase images, measures of 
velocity can be obtained while the signal phase changes caused by 
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susceptibility and field inhomogeneity are reduced or removed from the 
velocity image. By applying a bipolar gradient pulses, Gv, as taught in paper 
by Moran in J. Magnetic Resonance Imaging, 1982, a signal phase is 
proportional to velocity of spins moving along the applied gradient axis. 
This pulses Gv are labeled G s v in Fig. 13. The phase of the signal is also 
determined by the linear added effect of time spent in a local susceptibility 
gradient and in the inhomogeneity of the main magnetic field, together 
called phase errors, Ep. The Ep can be nulled by encoding the two images 
with equal but opposite velocity phase shift and combining the two images 
so that the velocity phases are coherently combined while the Ep are nulled. 
One approach is to apply a Gv pulse between the first and second RF 
excitations and a second Gv pulse after the second RF excitation where the 
first Gv pulse has twice the gradient amplitude and opposite polarity, -2Gv. 
The image from the first RF excitation has a net phase of (Gv -2Gv) equal 
to -Gv. The image created by the second RF excitation does not experience 
the first -2Gv pulse and has a net phase contribution by only the second 
gradient pulse, +Gv. The net phase shifts in the images also include the 
effects of Ep. A linear complex combination of the data from the two 
acquired images nulls the Ep contribution to phase and leaves a combined 
image with only velocity dependent phase (Pv), where (Imagel (-Pv +Ep) - 
Image2 (+Pv +Ep) = Image3 (Pv) . A final image has velocity phase (Pv) 
which is corrected for Ep and is dependent on Gv and the duration of Gv. 
Many different combinationsof Gv pulses on Gx , Gy and Gz axis can be 
applied, so that several final Ep corrected images are obtained. In fact, the 
velocity vector components Vx, Vy, and Vz can be obtained from the same 
echo train sequence, all of Which are corrected for Ep. This would require as 
many as 6 RF excitation pulses to create 6 images with different values of 
Pv dependent on Gv applied on three different gradient axis. An 
appropriate linear combination of the acquired images gives a phase image 
obtained from the complex data to display Pv as a measure of velocity. 
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Other encoding methods can be applied as gradient pulses between 
the RF excitation pulses. One example is illustrated in Fig. 14. 
Combinations of either Gc or Gv applied between RF pulses gives 
corresponding images of a spatial region obtained simultaneously from a 
5 single echo train. Also, the sequence can be repeated in iterative cycles (c) 
during different TR periods with varying either the Gp phase encoding pulse 
areas, or with different Gc or Gv pulses. The additional cycles of the 
sequence can be combined to make higher spatial resolution. Yet a different 
useful variant is for the additional c cycles of the sequence to be used to 

10 make additional images of the same slice region, each image using a 
different Gc or Gv so that additional correction images are obtained and 
used to be combined into the final image for improved correction of 
susceptibility or measurement of velocity. 

Conditioning gradient pulses can be used for purposed other than, or 

15 in addition to, reducing or eliminating susceptibility dropout. Fig. 1 illustrates 
an example that is otherwise similar to Fig. 3 but omits the conditioning 
gradient pulses of Fig. 3 and applies different phase-encoding to the same 
slice volume by means of a conditioning gradient pulse GpBlip applied on 
the Gp axis, in the time between the two RF pulses. An advantage of this 

20 variant compared to a conventional EPI pulse sequence lies in the time 
saving. If a conventional EPI sequence uses N phase-encoding lines (e.g. 
., N=128), then it has to switch the gradients at least N-1 times. The variant 
illustrated in Fig. 7, which uses two excitations, each differently phase- 
encoded, reduces the number of gradient switchings by more than a factor 

25 of 2, to be precise, to N/2-1 . 

Fig. 8 illustrates a similar approach but applied to 3D phase-encoding 
partitions. Here the gradient G3dP1 is placed on the slice-select axis Gs, 
between the two RF pulses, and an opposite polarity gradient G3dP2 is on 
the same Gs axis, but follows the second RF pulse. The sequence variant 

30 for 3D-encoding makes it possible to acquire two (or more) different 3D 
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Phase-encoding steps within one echo train, thus yielding a shorter 
measurement time. 

in .he axamp.es o. the new process discussed above. ,he echo sampling 
. surname, i.e. if the acquired K-space da.a is cu. in ha ,f, ,he echoes are fa,„ng 
> «o me m.ddle o, each par,. However. „ is possible .o use an asymme.ric 
readout, which allows separating the echoes within the k-space (moving t hem 
further apart). This reduces the interference of (he high frequency components of 
the echoes and can be used to minimize image artifacts. To achieve this 
separation the pre-phasing gradients in the read direction can be changed 
Therefore, an asymmetry-factor k is introduced, which is 0 for symmetric readout 
and 1 for completely asymmetric readout, where the echo lies at the edge of 
acquired k-space. Then, the total momentum of Grt P and 6r2P is calculated as 
follows: 

totalMomentum(GrlP) = totalMomentum(Greadout) *1/2 *(i +k ) 
totalMomentum(Gr2P) = totalMomentum(Greadout) *l/4 *(1-k). 

Fig. 9 illustrates echo spacing for three k-values, and shows that echoes 
are centered in k-space at k=0 but move further apart from each other as the 
value of k increases. 

As illustrated in Fig. 10, phase-cycling can be used to separate the signals 
of the RF pulses. This allows the echoes of all RF pulses to fall on top of each 
other. This corresponds to an asymmetry-factor k=-1 . But. multiple acquisitions 
are performed to phase-cycle the RF pulses, which places limits on the 
advantage of the method to acquire multiple images with different conditioning at 
the same time. The information can be acquired at the same time but averaged 
over the multiple acquisitions. In standard imaging techniques the information is 
not averaged over time, because separate readout periods are used, and the 
s-gnals are acquired to different time points, which can yield significant artifacts 
due to motion, physiological changes, etc. To increase time efficiency, the 
different phase cycles can be encoded within one echo train, by using both 
Phase-cycling and S.R-technology. As shown in Fig. 10. two RF pulses a1 and 
a2 are applied, which have different zshim-values, according to the gradient 



16 



WO 03/098233 



PCT/US03/15659 



pulse GsC between them. The difference here is, that no read gradient is applied 
between a1 and a2, thus making the echoes fall on top of each other and forming 
the signal s12. After those two RF pulses a1 and a2, a defocusing gradient pulse 
separates the combined signals s12 from the next two signals s34, which result 
5 from the RF pulses a3 and a4. The phase of a4 is changed be 180° compared to 
a2. This allows the process to separate the signals s1 and s3 from the signals s2 
and s4 by simply adding or subtracting the divided k-space data, respectively. 
Thus in turn allows the process to suppress signal leakage from one MR image 
to another image that has different conditioning. Differences in the phase 

10 evolution due to the slightly different echo times can be addressed by adding an 
additional pre-phasing gradient pulse GpP, which moves the echo of the RF 
pulses a1 and a2 to an earlier line in k-space than the echo of a3 and a4. For 
this, the distance between a1 and a3 (the same as the distance between a2 and 
a4 is a multiple of the time needed to acquire a single line in EPI readout. This 

1 5 factor gives the number of phase-encoding lines the echoes of a1 and 2 have to 
be shifted up. This method shows some similarity with POMP, where two 
different slices are acquired together and are separated by phase-cycling, too. 
However, unlike POMP, in this variant of the new process the summed signals 
originate from the same slice volume. 

20 The new process is not limited to the type and number of conditioning 

gradient pulses used in the examples discussed above. For example, instead of 
using only one conditioning gradient pulse Gc1 , which is located between the two 
excitation pulses RF1 and RF2, it is possible to add another conditioning gradient 
pulse Gc2 after RF2. In general, the magnetization that is excited by RF1 "feels" 

25 the sum of both conditioning gradients Gc1 and Gc2, while the magnetization 
that is excited by RF2 only "feels" the effect of Gc2. 

In some of the above examples two excitation pulses RF1 and RF2 are 
used to produce two independent images with different conditioning. However, 
the new process is not limited to only two simultaneously refocused images. One 

30 way is to add more RF excitation pulses to produce more primary echoes, which 
can be encoded according to the described method for two excitation pulses 
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(addition of pre-phasing gradient pulse along the read-direction, called Gr1P 
above). However, the signal-to-noise ratio is likely to be reduced since more RF 
pulses are applied on the same slice volume. An ISMRM abstract that the 
named invention have prepared for publication, but which is not published as of 
the filing date of this patent application, describes how three images can be 
encoded with only two excitation pulses RF1 and RF2. The second RF pulse 
RF2 will not only produce a primary echo, but also act as a refocusing SE-echo 
pulse on the magnetization of the first RF pulse RF1. This will yield a spin echo 
after RF2, which can be used for image acquisition, too. Therefore, not only two 
but three images can be acquired, which have different conditioning, such as 
different z-shim correction values. However, the amount of conditioning of the 
three images is not independent any more. Since the magnetization, which forms 
the spin echo, was reversed by RF2, the phase history of those spins is mirrored, 
i.e. using two conditioning gradient pulses Get and Gc2 after application of RF1 
and RF2, respectively, the primary echo of RF1 experiences the effect of 
Gc1+Gc2, while the primary echo of RF2 is prepared with Gc2 only. However, 
the spin echo of RF1 , which is produced by RF2, "feels" -Gc1+Gc2. If the values 
of Gc1 and Gc2 are adjusted appropriately, three differently conditioned images 
can be acquired with two RF excitation pulses. The inventioned can be extended 
to use multiple RF excitation pulses, more than two pulses, to create pleural 
signals. The invention can be combined with additional pre conditioning RF 
pulses and gradient pulses which eliminate lipid signal and leave only water 
signal in the image. 

A novel EPI sequence described in the ISMRM abstract that is yet 
unpublished is used to eliminate susceptibility dependent signal drop out in the 
brain, limiting to fMRI. The process simultaneously acquires 3 different gradient 
compensated images at the identical slice plane or volume and essentially at the 
same time. The net combined single shot T2* SIR EPI image recovers the signal 
lost due to susceptibility in the temporal and frontal lobes of the brain. 

Echo planar imaging is extensively used for functional MRI and brain 
activation mapping. EPI has well characterized susceptibility artifacts causing 
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signal drop out in brain regions adjacent to air tissue interfaces in the inferior 
temporal lobes, gyrus rectus and orbitofrontal lobes. This problem of incomplete 
brain coverage is perhaps the major imaging problem remaining in brain 
activation studies. Several methods have been developed to recover 
5 susceptibility signal loss by means of applying magnetic gradients to correct for 
local field inhomogeneity. These methods require repeating the EPI sequence 
with multiple excitations with different gradient correction pulses, or z-shim. This 
lengthens the imaging experiment and imposes further constraints in performing 
single trial studies. A novel EPI sequence incorporates the Simultaneous Image 

10 Refocusing (SIR) technique to overcome susceptibility artifacts in single shot EPI 
while maintaining T2* weighting for BOLD contrast. The earlier described SIR 
technique acquires multi-slice images in a single echo train using multiple 
excitation RF pulses, simultaneously refocusing all signals with each polarity 
switching of the read gradient. It is modified to simultaneously obtain multiple 

1 5 images of the same slice volume in one echo train. The simultaneously acquired 
images are encoded with different gradient correction pulses and combined to 
form a single shot EPI image with recovery of the susceptibility signal drop out in 
the brain. 

Fig. 11 illustrates this sequence. Two slice selective RF pulses are applied 
20 to the same slice location. Three signal components are generated as the result, 
namely 2 FIDs, one from each RF pulse and an echo generated by the two RF 
pulse. These signals share the same phase encoding, but are separated along 
the readout direction. Each signal component receives a different net area of the 
gradient pulses for susceptibility correction, applied on the slice direction (Gs). 
25 With a single shot data acquisition, three complete images are obtained. These 
images are combined after normal 2D FT image reconstruction to make an 
image without signal dropouts. 

The sequence was implemented on a Siemens 1.5T Sonata scanner. 
Typical imaging parameters were: single shot, T2* TE = 50 msec, 256 x 64 
30 matrix, which is separated into three 85 x 64 images, FOV = 30 cm, slice 
thickness = 5 mm. Fig. 12 shows in-vivo images of a healthy volunteer chosen 
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from a multi-slice acquisition in regions of most severe signal drop-out. The top 
row shows images reconstructed using only signals from FID 1, which does not 
receive any correction gradient, essentially a conventional T2* EPI. The Lower 
row shows the SIR EPI images using all three components combined, including 
5 FID 1. In comparing the signal drop out regions, it is critical to note that the EPI 
(top row) and SIR EPI (lower row) have identical T2* weighting at TE of 50msec. 

The multi-slice SIR EPI sequence successfully reduces or eliminates 
susceptibility dependent signal drop out in the brain's orbito-frontal cortex and 
temporal lobes, which limits conventional EPI sequences with the same T2* 

10 weighting. SIR EPI, therefore, provides greater image coverage of the brain as 
needed for many fMRI experiments. It is significant that the signal drop out 
regions are recovered in a single shot sequence, unencumbered by longer 
acquisition times from sequence repetitions needed by other susceptibility 
correction schemes. The new single shot SIR EPI sequence provides intrinsic 

15 gradient correction and recovery of signal in brain regions of high susceptibility, 
limiting EPI sequences. SIR EPI therefore may likely provide significant 
improvements over conventional EPI in functional MRI experiments. 

Hard (non-selective) RF pulses can be used to reduce the time spent for 
excitation and thus reduce the time distance (and phase evolution) between RF 

20 pulses. This is similar to BURST-like techniques, where multiple hard RF pulses 
are used with different phase encoding steps to form an image. In contrast, the 
new process can use the hard RF pulses not to form one image but to form 
multiple images with different conditioning, which is performed between (and 
after) the RF excitation pulses. To achieve slice selectivity, saturation pulses can 

25 be applied, to suppress ad but the wanted slice, in a manner similar to that 
proposed in the Heid patent incorporated by reference herein. 

Persons skilled in this technology will understand that the above examples 
are not limiting to the scope of the disclosed process, and that many variations 
are possible within the teachings set forth above, and that MRI systems 

30 configured to carry out the disclosed processes, such as by suitable pulse 
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sequence or software changes, constitute new machines that are within the 
scope of the disclosure herein. 
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What is claimed: 

1. An MRI method comprising: 
subjecting a body to a main magnetic field; 

applying to the body a first RF excitation pulse and a first magnetic field 
gradient to define a vollume in the body, and later applying to the body a second 
RF excitation pulse and a second magnetic field gradient to define the same 
volume; 

giving the volume defined by the first RF pulse and first gradient a phase 
history different from that of the same volume defined by the second RF pulse 
and second gradient 

obtaining first MRI signals for the volume as excited by the first RF pulse, 
and second MRI signals for the same volume as excited by the second RF pulse 

said first and second MRI signals alternating in time; 

reconstructing a first MRI image from the first MRI signals, and a second 
MRI image for the same slice from the second MRI signals; 

combining said first and second MRI images into a net image; and 

using said net image. 

2. An MRI method as in claim 1 , in which said different phase histories 
are selected to give one of the first and second images a correction of 
susceptibility different from any such correction of the other image. 

3. A method as in claim 1 in which said gradients are slice selecting 
gradients 

4. An MRI method comprising: 

obtaining substantially simultaneously MRI signals for two images of the 
same volume encoded with a different weighting or rephasing by correction 
magnetic gradient pulses; 
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combining information derived from said MRI signals into an image that 
has signal recovery in regions of high susceptibility phase errors. 
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